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Abstract 

Background: Transforming growth factor-beta 1(TGF-[S1) is involved in the development of acute rejection (AR) episodes in 
solid organ transplant recipients; and a number of studies have been conducted to investigate the combined effects of 
human TGF-pi gene (TGFB1) +869 T/C and +915 G/C polymorphisms on AR risk. However, the results obtained are 
inconclusive. 

Methods: Eligible studies that investigated the haplotypic association between TGFB1 +869 T/C and +915 G/C 
polymorphisms and AR risk were comprehensively searched in the PUBMED, EMBASE, China National Knowledge 
Infrastructure, and Wanfang Database. Statistical analyses were performed by using STATA 12.0 and Review Manager 5.0. 

Results: Fourteen eligible studies with 565 AR cases and 1219 non-AR cases were included. Overall, a significantly decreased 
risk was detected in patients carried with intermediate producer (IP) haplotypes (T/C G/C, T/T G/C, and C/C G/G) and/or low 
producer (LP) haplotypes (C/C G/C, C/C C/C, T/T C/C, and T/C C/C) compared with high producer (HP) haplotypes (T/T G/G 
and T/C G/G; IP vs. HP: OR = 0.75, 95% CI, 0.58-0.96, P heterogeneity =0.238; IP/LP vs. HP: OR =0.77, 95% CI, 0.61-0.98, 
P heterogeneity =0.144). In addition, subgroup analysis by transplant types demonstrated a similar association in patients 
receiving heart transplant (IP vs. HP: OR =0.32, 95% CI, 0.14-0.73, P heterogeneity =0.790; IP/LP vs. HP: OR =0.41, 95% CI, 

0.20-0.85, P heterogeneity =0.320). 

Conclusions: The current meta-analysis and systematic review indicated that recipient TCFB1 HP haplotypes were 
significantly associated with an increased risk for AR in solid organ transplant recipients, particularly patients receiving 
cardiac allograft. 
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Introduction 

Transforming growth factor-beta 1 (TGF-pi) is a multifunc- 
tional cytokine ubiquitously produced by a wide variety of cells, 
including T lymphocytes, monocytes, vascular endothelium, and 
renal tubular cells [1]. Functionally, TGF-fH has been proven to 
be of fundamental importance in the development of various 
disorders [2] , including coronary heart disease [3] , human cancers 
[4], rheumatoid arthritis [5], asthma [6] and transplant rejection 
[7,8]. In the setting of solid organ transplants, TGF-fil has been 
conventionally recognized as a guardian against acute rejection 



(AR), as higher level of TGF-(3l in the graft tissue and serum was 
found in non-AR recipients than those suffering AR [9—11]. 
However, several novel lines of evidence have challenged the 
beneficial effects of TGF-fil on transplant recipients [12,13]. 
Although the functional role of TGF-|3l in AR initiation remains 
elusive, this cytokine is believed to exert pivotal and complicated 
functions in AR episodes. 

The human TGF-|3l gene (TGFB1) is mapped on the 
chromosome 19ql3. 1-13.3 with seven exons and six introns, 
whose regulation and expression is influenced by the presence of 
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single nucleotide polymorphisms (SNPs) [14]. Among these SNPs, 
+869 T/C (also known as rs 1800470, T29C, or LeulOPro) and + 
915 G/C (also termed as rs 1800471, G74C, or Arg25Pro) 
polymorphisms in the first exon of TGFB1 have been the focus 
of extensive researches and donor TGFB1 +869 T/C polymor- 
phism has been proven to be significantly associated with AR risk 
[15]. These two SNPs together result in nine potential inherited 
haplotypes, which could be categorized into three groups 
according to the in vitro production levels: high producer (HP) 
(T/T G/G and T/C G/G), intermediate producer (IP) (T/C G/ 
C, T/T G/C, and C/C G/G) and low producer (LP) (C/C G/C, 
C/C C/C, T/T C/C and T/C C/C) [16,17]. Since the first study 
conducted by Pelletier et al to evaluate the haplotypic association 
of TGFB1 +869 T/C and +915 G/C polymorphisms with AR risk 
in kidney transplant recipients [18], numerous molecular epide- 
miological studies have been conducted in different solid organ 
transplants, including kidney transplants [19-26], liver transplants 
[9,27,28], and heart transplants [13,29]. However, the results of 
these studies were inconclusive. 

In this meta-analysis, we integrated the data from all eligible 
studies to explore 1) the combined effects of recipient TGFB1 +869 
T/ C and +9 1 5 G/ C polymorphisms on AR risk after solid organ 
transplantation and 2) the potential influence of covariants such as 
ethnicity, transplantation types, and immunosuppressive protocols. 

Materials and Methods 

Identification of eligible studies 

This meta-analysis was conducted and reported in accordance 
with the PRISMA (Preferred Reporting Items for Systematic 
reviews and Meta- Analyses) guidelines (Checklist SI) [30]. To 
identify all eligible studies that investigated the haplotypic 
association of TGFB1 +869 T/C and +915 G/C polymorphisms 
with AR risk in solid organ transplantation, a comprehensive 
electronic search of PUBMED, EMBASE, China National 
Knowledge Infrastructure (CNKI), and Wanfang databases was 
performed until November 29, 2013. To search and include as 
many related studies as possible, we applied various combinations 
of the following medical subject headings and key words: 
transforming growth factor beta-1, TGFfS-1, or TGFB1; acute 
rejection, early allograft outcome, or graft rejection; transplanta- 
tion or organ transplantation; polymorphisms or variants. 
Furthermore, the reference lists of reviews and retrieved articles 
were manually screened for additional studies No restrictions were 
placed on language, and only published studies with full-text 
articles were included. 

Inclusion and exclusion criteria 

The studies identified from the above mentioned databases were 
screened by two independent authors (Yu-Zheng Ge and Ran Wu) 
according to the following predesigned inclusion criteria: 1) case- 
control design; 2) evaluating the correlation of TGFB1 haplotypes 
(at position +869 and +915) with AR risk; and 3) providing 
sufficient data to calculate the odds ratio (OR) and its 
corresponding 95% confidence interval (CI). When several studies 
with overlapping data were eligible, those with smaller sample size 
or less reliability were excluded. In addition, studies without 
detailed information were excluded, after the efforts to extract data 
from the original paper or contact the corresponding authors 
failed. 

Data extraction 

Data of eligible studies were extracted by two reviewers (Yu- 
Zheng Ge and Tian-Ze Lu) independendy and in duplicate 



according to the predesigned data-collection form. The following 
information was extracted: last name of first author, publication 
year, country of origin, ethnicity, transplantation type, immuno- 
suppressive protocol, number of both ARs and non-ARs and 
phenotypic distribution in both groups. Different ethnic descents 
were categorized as Asian, Caucasian, African and Mixed (which 
included more than one ethnic descent). Transplantation types 
were characterized as renal, liver, and heart transplantation. Based 
on different calcineurin inhibitors (CNIs) within the immunosup- 
pressive protocol, the studies were divided into two subgroups: 
Cyclosporine A (CsA) group and CsA/tacrolimus (FK506) group. 
Discrepancies occurring during the process of study inclusion and 
data extraction were resolved by discussion with a third reviewer 
(Wen-Cheng Li), and consensus on each item was achieved 
eventually. 

Statistical analysis 

Crude OR and its corresponding 95% CI were used to assess 
the strength of association between the haplotypes of +869 T/ C 
and +915 G/C polymorphisms in TGFB1 gene and AR risk. The 
HP haplotypes (T/T G/G and T/C G/G) were used as the 
baseline for calculation of ORs in three different comparisons (LP 
vs. HP, IP vs. HP and LP/IP vs. HP). Stratified analyses were also 
conducted based on ethnicity, transplant types, and immunosup- 
pressive protocols. The statistical significance of the pooled OR 
was assessed with the Z test, and P<0.05 was considered 
significant. 

Chi-square based Q;test was performed to measure between- 
study heterogeneity, and the presence of heterogeneity was 
considered significant if P<0.10 [31]. When the heterogeneity 
was absent, the fixed-effect model (Mantel-Haenszel method) was 
used to pool the data from different studies [32]; otherwise, the 
random-effects model (DerSimonian and Laird method) was 
applied [33]. Sensitivity analyses were conducted to identify the 
effect of individual study on pooled results and to test the reliability 
of results by deleting a single study each time [34] . To determine 
the presence of publication bias, both Begg's funnel plot and 
Egger's linear regression test were conducted, and P<0.05 was 
considered significant [35,36]. 

All statistical analyses were performed with STATA software 
(version 12.0; Stata Corporation, College Station, Texas, USA) 
and Review Manager (version 5.0; Cochrane Collaboration, 
Oxford, UK). 

Results 

Characteristics of eligible studies 

Fourteen eligible articles were identified according to the 
predesigned selection criteria, and the detailed screening process 
was shown in Figure 1. Among the 14 eligible studies, 9 were renal 
transplantation [18-2 6], 3 were liver transplantation [9,27,28], and 
2 were heart transplantation [13, 2 9]. As for ethnicity, 7 were 
studies of Caucasians [9,22,23,27,28,37], 3 studies were of Asians 
[19-21], and 4 studies of mixed ethnicity [13,18,26,29]. Regarding 
the immunosuppressive protocols, patients were uniformly pre- 
scribed with CsA in 5 studies[18-20,23,29] while in the remaining 
9 studies were applied with either CsA or FK506 [9, 13,2 1,22,24- 
28] (Table 1). 

Quantitative data synthesis 

A total of 14 studies including 565 AR and 1219 non-AR cases 
were identified to assess the haplotypic association of TGFB1 +869 
T/ C and +9 1 5 G/ C polymorphisms with AR risk of solid organ 
transplantation. Overall, a significantly decreased risk was 
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detected in two comparisons (IP vs. HP: OR =0.75, 95% CI, 
0.58-0.96, P heterogeneity = 0.238, Figure 2; IP/LP vs. HP: OR 
= 0.77, 95% CI, 0.61-0.98, P heterogeneity =0.144, Figure 3). 
Subgroup analyses based on ethnicity, transplant types, and 
immunosuppressive protocols were subsequendy conducted, and 
the results demonstrated a remarkably decreased risk in heart 
transplant recipients (IP vs. HP: OR =0.32, 95% CI, 0.14-0.73, 
P heterogeneity =0.790, Figure 2; IP/LP vs. HP: OR =0.41, 95% 
CI, 0.20-0.85, P heterogeneity = 0.320, Figure 3) and in patients 
uniformly administrated with CsA (IP vs. HP: OR = 0.57, 95% CI, 
0.34—0.94, P heterogeneity = 0.491). Table 2 represents the strength 
of association between TGFB1 haplotypes (at position +869 and + 
915) and AR risk in transplant recipients. 

Heterogeneity test and sensitivity analysis 

No significant between-study heterogeneity was observed in this 
meta-analysis (Table 2). Sensitivity analysis was also performed to 
explore the potential influence of each individual study on the 
overall results by deleting one single study each time from the 
pooled analysis. NO substantial change was demonstrated in the 
overall studies, indicating that no individual study could affect the 
pooled OR significantly (data not shown). 



Publication bias 

To examine the publication bias of the currently available 
literature, both Begg's funnel plot and Egger's test were 
conducted. The shape of the funnel plots did not reveal any 
evidence of obvious asymmetry in all comparison models (Figures 
SI, S2, S3). Then, the Egger's test was used to provide statistical 
evidence for funnel plot symmetry. The results also did not show 
any evidence of publication bias (Table SI). 

Discussion 

Organ transplantation has been recommended as the most 
optimal treatment choice for patients suffering end-stage disease. 
However, AR still remains a crucial determining factor which 
influences the short-term function and long-term outcome of both 
recipients and allografts [38-40] . To the best of our knowledge, 
this is the first meta-analysis focusing on the combined effects of 
TGFB1 +869 T/C and +915 G/C polymorphisms on AR risk. In 
the current study, we provide evidence that patients carried with 
TGFB1 HP haplotypes (T/T G/G and T/C G/G) are more likely 
to suffer from AR after solid organ transplantation (specifically 
heart transplantation), which could be utilized to identify patients 
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Figure 2. Forest plot for risk of acute rejection associated with TGFB1 haplotypes (IP vs. HP) stratified by transplant types. For each 
study, the estimate of OR and its 95% CI is plotted with a box and a horizontal line. Filled diamond pooled OR and its 95% CI. 
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predisposed to AR and potentially benefiting from tailored 
immunosuppressive protocol. 

After pooling the data from 14 eligible studies, we shown that 
patients with TGFB1 IP and/ or LP haplotypes were less likely to 
suffer from AR than those with HP haplotypes (IP vs. HP: OR 
= 0.75, 95% CI, 0.58-0.96; IP/LP vs. HP: OR =0.77, 95% CI, 
0.61-0.98). Considering that ethnic background could influence 
the frequency of genotypes and haplotypes in some context, we 
conducted subgroup analysis based on ethnicities but failed to 
detect any significant association. Further stratified analysis by 
transplant types demonstrated a remarkably decreased AR risk for 
heart transplant recipients carried with TGFB1 IP and/or LP 
haplotypes, which was in agreement with the study conducted in a 
cohort of 1 1 1 pediatric cardiac transplant recipients [13]. With 
respect to various CNIs applied in suppressive protocols, CsA and 
FK506 have been proven to influence the serum level of TGF-(3l 
differentially [41,42]; thus, we intended to investigate the various 
impacts of different CNIs on the relationship between the TGFB1 
haplotypes and AR risk. However, FK506 has never been 



administrated as a unique CNI in any included studies. Therefore, 
we divided the studies into two subgroups: CsA arm and CsA/ 
FK506 arm, and then conducted a subgroup analysis. As a result, 
patients with TGFB1 IP haplotypes were less likely to be affected 
by AR in CsA arm (IP vs. HP: OR =0.57, 95% CI, 0.34-0.94), 
which could be interpreted as a clue that the recipients with 
TGFB1 IP haplotypes can benefit from CsA after organ 
transplantation. 

TGF-fH is a pleiotropic and multifunctional cytokine with 
immunosuppressive and fibrogenic properties, which may play a 
central role in both the initiation and propagation of AR and 
chronic rejection (CR) [7,8]. The expression of TGF-fSl has been 
proven to be regulated by the two SNPs (+869 T/C and +915 G/ 
C polymorphisms), whose combination could be divided into three 
groups (LP, IP, and HP). However, the results obtained have 
challenged the conventional concept that TGF-fSl may inhibit the 
initiation of AR episodes. Our team previously demonstrated that 
donor TGFB1 +869 T/C HP genotype (TT) was significandy 
associated with decreased AR risk [15]. Considering the obvious 
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Figure 3. Forest plot for risk of acute rejection associated with TGFB1 haplotypes (LP/IP vs. HP) stratified by transplant types. For 

each study, the estimate of OR and its 95% CI is plotted with a box and a horizontal line. Filled diamond pooled OR and its 95% CI. 
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discrepancies, several questions have been raised: 1) The impact of 
various haplotypes of TGFB1 +869 T/C and +915 G/C 
polymorphisms on actual production levels of TGF-J31. The 
effects of the nine potential haplotypes on TGF-pi expression 
levels have been previously studied in vitro and among patients 
without organ transplant [16,43,44]. However, no significant 
association between TGFB1 SNPs and TGF-pi plasma levels or 
intragraft mRNA levels was detected in a cohort of renal 
transplant recipients [45] . 2) The influence of various immuno- 
suppressive protocols on the expression level of TGF-f11. CsA and 
FK506 were found to influence the serum level of TGF-pi 
differentially [41,42,45], while the impact of other immunosup- 
pressive regimens still remains elusive. 3) The potential effects of 
donor and/ or donor-recipient pair TGFB1 haplotypes on AR risk. 
The current study is confined to recipient SNPs; however, in some 
context, donor polymorphisms may contribute much more to AR 
risk than recipient SNPs [15,46]. 4) The exact etiology of AR 
episodes and the biological functions of TGF-pi in the develop- 
ment of AR still remain unclear. 



The current meta-analysis focused on the combined effects of 
TGFB1 +869 T/C and +915 G/C polymorphisms rather than one 
single SNP on AR risk [15], which could help to derive a precise 
estimation of the roles of TGFB1 SNPs in the development of AR 
episodes. However, several limitations should be considered when 
interpreting the results. First, considering the unavailability of 
other detailed information, we did not conduct stratified analysis 
based on some cofactors such as follow-up time, gender, age, panel 
reactive antibodies level, human leukocyte antigens mismatch and 
donor source, which may influence the results. Second, the limited 
number of AR cases and non-AR cases may lead to a relatively 
small power. Third, only published studies with sufficient data 
were included, thus, publication bias may have occurred even 
though results of both Begg's test and Egger's test did not detect it. 
Last but not least, the meta-analysis is retrospective due to the 
methodological limitations. 

In summary, this meta-analysis suggested that recipient TGFB1 
HP haplotypes of +869 T/C and +915 G/C polymorphisms (T/T 
G/G and T/C G/G) might be a possible genetic susceptibility 
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locus for AR after solid organ transplantation, which could be 
utilized to identify patients predisposed to AR and potentially 
benefiting from personalized immunosuppressive protocol. In 
addition, monitoring TGFB1 could help manage GR to some 
extent, as TGFB1 triggers fibrogenesis linked to chronic rejection 
(CR). Further well-designed and unbiased studies with larger sample 
size, diverse ethnicities, donor-recipient pairing and various 
applications of CNIs should be conducted to verify our findings. 
Furthermore, functional studies of TGFB1 gene polymorphism are 
warranted to understand the underlying mechanisms. 

Supporting Information 

Figure SI Begg's funnel plot for publication bias test 
(LP vs. HP for TGFB1 haplotypes). 

(TIF) 

Figure S2 Begg's funnel plot for publication bias test (IP 
vs. HP for TGFB1 haplotypes). 

(TIF) 
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